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Abstract 

Electrons and electron neutrinos in the inner core of the core-collapse supernova are 
highly degenerate and therefore numerous during a few seconds of explosion. In contrast, 
leptons of other flavors are non-degenerate and therefore relatively scarce. This is due 
to lepton flavor conservation. If this conservation law is broken by some non-standard 
interactions, electron neutrinos are converted to muon and tau-neutrinos, and electrons - to 
muons. This affects the supernova dynamics and the supernova neutrino signal in several 
ways. In particular, the total neutrino luminosity in the first second of the collapse is 
increased due to the larger free path of the non-electron neutrinos. This effect may have 
important consequences as the increase of the neutrino luminosity is known to facilitate 
the explosion. We consider an extension of the Standard Model by scalar bileptons (i.e. 
. heavy scalars with lepton number 2) which mediate lepton flavor violation. It is shown 

JZL|' that in case of TeV-mass bileptons the electron fermi gas is equilibrated with non-electron 

species inside the inner supernova core at a time-scale ~ (1 — 100) ms. In particular, a 
' scalar triplet which generates neutrino masses through the see-saw type II mechanism is 

considered. Non-observation of rare decays together with known data on neutrino mixing 
and neutrino masses restrict possible lepton flavor violation effects in this case. However a 
region in the parameter space of the model exists which fits all the experimental constraints 
' and provides lepton flavor violation sufficient for observable effects in supernova. 

m 
oo 

°^ ; 1 Introduction 

The present contribution is largely based on our recent paper [1]. In addition, here we present 
an extended discussion of the increase of the total neutrino luminosity during the first second 
of the collapse and possible consequences of this important effect for supernova dynamics (see 
section I5.3P . 

• Jh 
V- 

If the Standard Model (SM) is extended by a triplet of heavy scalar fields interacting with 
leptons and bosons of the SM, it becomes possible to generate Majorana neutrino masses through 
the see-saw type II mechanism [2]-[3]- Originally the mass of this triplet was assumed to be 
at a GUT scale; recently, however, TeV-scale see-saw became popular. In particular, it was 
argued that TeV-scale see-saw type II avoids electroweak hierarchy problem, encountered in 
GUT-scale see-saw, and meets the naturalness argument (see e.g. |5j|6j). See-saw II signatures 
at the LHC (see e.g. |7J |8j and references therein) and in low-energy experiments (see e.g. [6j [9j 
and references therein) are extensively investigated. 

In the present paper we point out that TeV-scale see-saw II affects the core-collapse super- 
nova physics, namely the supernova neutrino signal and the dynamics of the explosion. This 
conclusion is valid in a certain area of the parameter space of the neutrino mixing matrix. The 
effect is due to the lepton flavor violating (LFV) four-fermion interactions, which are induced 
by the scalar exchange. 
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Scalar triplet, which enters the see-saw II, is a particular case of scalar bilepton, i.e. scalar 
which couples to a pair of leptons. Another type of bilepton is a doubly charged singlet scalar. 
If possessing TeV-scale mass, it may also be produced at LHC and influence the results of low- 
energy experiments. We include the analysis of its effect on the LFV in supernova in the present 
paper, because it mediates LFV interactions analogously to the doubly charged component of 
the scalar triplet. 

It is known that LFV processes may alter drastically the conditions in the interiors of a 
core-collapse supernova [10] [H]. We discuss several possible consequences of such alteration, 
the most intriguing one being the possibility to facilitate the explosion. 

The paper is organized as follows. In section 2 we estimate the value of the effective four- 
fermion LFV coupling necessary to significantly affect the supernova physics. In section 3 we 
consider a simple extension of the SM with a doubly charged bilepton. Section 4 is devoted to 
see-saw type II extension of the SM. Influence of LFV inside the supernova core on supernova 
physics is discussed in section 5. The results are summarized in section 6. 

2 Model independent considerations 

An explosion of a non-thermonuclear supernova is a manifestations of a collapse of a progenitor 
star core. Although numerous studies have failed to reproduce an explosion of a core-collapse 
supernova, there exists a commonly accepted general picture of the collapse (which is referred to 
as "standard picture" in what follows; for the details see, for example, |12j). When the mass of 
the iron core of a massive star reaches the Chandrasekhar limit, the infall phase of the collapse 
starts. The core contracts due to the gravitational attraction. Some fraction of electrons is 
converted to electron neutrinos through the inverse beta processes. When the density of the 
inner part of the core reaches the nuclear density value, ~ 3 • 10 14 g/cm 3 , the infalling matter of 
the outer core bounces from it. A shock wave emerges, which starts to propagate outwards. A 
powerful shock could expel the stellar envelope and produce the supernova explosion. However, 
according to the detailed numerical calculations [13] , the shock wave looses its energy and stalls 
after some tens of milliseconds. There exist several scenarios of presumably successful explosion 
(e.g. [H]), but neither of them is justified to date by detailed self-consistent calculations with 
neutrino transport. 

During the collapse extreme values of density, temperature, electron and electron neutrino 
chemical potentials are reached. For the estimates we use the reference values of the supernova 
core parameters as presented in Table [TJ 
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Table 1: Typical conditions in the inner supernova core (m < 0.5M@) during the first 50 ms 
after core bounce as obtained by us with the help of the open-code program Boom [15]. Only 
SM interactions are taken into account. 

This data is obtained by us with the help of the open-code programme [15] (in which only 
SM interactions are taken into account), being in good agreement with the results of more 
sophisticated calculations, see e.g. [E][T7]. The presented values are the mean values for the 
inner supernova core (m < O.5M0); the central values are even greater. Note that in the stan- 
dard picture of the collapse electrons and electron neutrinos are plentiful and highly degenerate 
during the first hundreds of milliseconds. Their total number roughly equals the total number 
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of electrons in the core before collapse, and decreases slowly due to the leakage of neutrinos 
from the core. By contrast, there are very few non-electron neutrinos and muons. Non-electron 
neutrinos have negligible chemical potentials; they are created in pairs with their antiparticles, 
concentrations being proportional to T 3 . Muons have chemical potential approximately equal 
to [i e — fj, Ue , their concentration has an additional suppression factor e~ i < m ^~i L i^l T . Lepton flavor 
conservation law prevents electrons and electron neutrinos to turn into other types of leptons. 

This picture is modified if there exist LFV processes beyond the SM [K)J[llJ. If electrons 
and electron neutrinos are converted into leptons of other flavors inside the inner cor 60, the 
dynamics of the collapse and the supernova neutrino signal is different from the SM case. In 
particular, the energy possessed by leptons is redistributed between all neutrino flavors, and 
degenerate seas of i>a and v T emerge. The following LFV reactions could be in principle relevant 
for supernova!! 
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The reactions are divided into two groups. In the first group all individual lepton flavor numbers, 
L e ,L^ and L T are changed by two units or unchanged at all (e.g. for the first reaction AL e = 
—2, AL M = 2, AL T = 0). In the second group at least some of the individual lepton flavor 
numbers are changed by one unit. As is discussed below, processes from the second group 
are tightly constrained by non-observation of rare decays and hardly can affect the supernova 
physics; in contrast, the processes from the first group can affect the supernova physics and at 
the same time fit the present experimental constraints. 

Let us estimate the value of the cross sections of the LFV reactions potentially important 
for supernova core. For this purpose we demand that such reactions should lead to the thermal 
equilibrium between e and v e on the one hand and leptons of other flavor (-s) on the other hand 
at a certain timescale At (note that equilibration between charged leptons and neutrinos of the 
same flavor is established virtually instantly due to the ordinary weak interactions). We consider 
two reference values for At. The first one, At\ = 1 ms, is characteristic for the dynamics of the 
supernova (e.g. the shock propagates through the core during several ms). The second one, 
At2, should characterize the Kelvin- Helmholtz timescale of the cooling of a proto-neutron star 
(it is somewhat larger than the neutrino diffusion timescale, see section UTT]) . which determines 
the long-term supernova neutrino signal. Although the inner core is deleptonized due to the 
diffusion at a timescale ~ 5 s [18] , we conservatively chose At2 = 300 ms, as at greater times 
the numbers presented in Table 1 are significantly changed. 

From the above-mentioned condition we find for the processes in which electrons are present 
in the initial state 

a > \/{Atn e c) ~ 3 • 10~ 48 (300 ms/At)cm 2 . (3) 

If only neutrinos are involved in the process, the cross section should be several times greater, 
as the neutrino density is several times smaller. The estimated value appears to be rather small 



Not to be confused with ordinary resonant neutrino flavor conversion due to the MSW-effect, which take 
place in the envelope. It affects the neutrino signal but do not affect the supernova dynamics. 
2 We consider processes in which only leptons take part. 
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compared to ordinary weak cross sections. For example, 

a{v e n -»• e~p) ~ 2 • 10~ 39 cm 2 (4) 

for 160 MeV neutrinos [13] . 

Reactions ([1]) and ([2]) at low energies are described by effective four-fermion terms in the 
Lagrangian, the cross section being of order of G^py^ 2 ^, where H e ,v e stands for the chemical 
potential of electron and electron neutrino correspondingly (every reaction is characterized, in 
general, by its own Glfv)- From ([3]) one gets a rough estimate for Glfv necessary for large 
LFV in supernova. In case of electron-electron collisions it reads 

Glfv > 4 • 10~ 4 v / 300 ms/Ai TeV~ 2 . (5) 

In case of neutrino-electron and neutrino-neutrino collisions the bound is somewhat larger due 
to the smaller neutrino density and chemical potential. Such bounds on Glfv ~ A 2 /M 2 indicate 
the mass M of intermediate heavy boson of order of several TeVs if the dimensionless coupling 
constant A is of order of 0.1. 

Are such values of Glfv compatible with present experimental constraints? This question 
is accurately treated in the following two sections, here we provide a preliminary discussion. 
The answer is different for two groups of reactions, (pQ) and ((21). In the former case, the best 
experimental limit stems from the non-observation of oscillations of muonium atom, /i~e + 
/x + e _ [2Q], and reads Glfv ^5 10 _1 TeV~ 2 . Although muonium oscillations involve only charged 
leptons, due to the SU(2) symmetry the bound on the effective constant is valid for the first 
three processes in ([T]), including those which involve neutrinos. As for the fourth process in ((TJ), 
it evades even this not very restrictive bound. Thus there is a room for LFV in supernova due 
to four reactions ((TJ) with |AL ejA1]T | = 0,2. 

As for reactions with \ALi\ = 1 d2]), they are tightly constrained by the rare lepton decays: 

fM — > eee, 

H -»■ 67, (6) 
r -»■ III. 

The constraints span from Glfv < 10" 2 TeV~ 2 (LFV decays of r) to Glfv < 10~ 5 TeV" 2 
(/i~ — > e~e~e + ), see [21 j and references therein. Therefore processes ([2]) can fit condition 

only marginally or can not fit it at all. Thus they are likely to be irrelevant for LFV in 
supernova. 

To conclude this section, we are interested in such an extension of the SM, in which four- 
fermion LFV processes with |ALj| = 1 are suppressed compared to the processes with |ALj| = 2. 
In the next section we give a simple but instructive example of such extension. In section 4 we 
show that see-saw type II model fits this requirements for certain values of neutrino masses and 
mixing parameters. 

3 Doubly charged scalar singlet 

Consider the SM extended by one doubly charged scalar A, singlet under SU(2). It can not 
be coupled to quarks due to its electric charge. The only possible coupling of such scalar to 
leptons reads 

c A = j2 x wWR A + h - c - ( 7 ) 

1,1' 

Here 

1 and I' are charged leptons, 1,1' = e,fj,,r; 

A/;/ is a matrix of coupling constants, summation over lepton flavors is implied. 
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Figure 1: ee — )■ fi/i LFV transition mediated by the doubly charged bilepton A 

C-conjugated right-handed fermions here and in what follows are defined as i^ c R = {ipR) c = 
7 2 (V'-R)*) standard representation of gamma-matrices being used. 

Let us choose the coupling matrix to be proportional to unit matrix: A//' = X5w. In this 
model low-energy processes with |ALj| = 1 are forbidden H 

The exchange of doubly charged bilepton gives rise only to one LFV reaction, interesting in 
the context of supernova core conditions (see Fig. [1]): 

e~e~ — > piT [f . (8) 

The corresponding effective low-energy LFV term reads 

^2 

C LFV = JjT e R e R^RVR)\ (9) 

which gives (see e.g. [22] [23] [24] ^ 

A 4 

cr(ee -)• fi/i) = v^(l - 2m 2 / s)\t max - t min \/8n. (10) 

Note that one should take into account the combinatoric factors and the identity of fermions 
in the final state while calculating this cross section. In the center-of-mass reference frame 

\tmax ~ tmin\ = ^E ey J ' E\ - m 2 . Taking E e ~ \i e ~ 2m ^ one obtains 

^^^" 2 - 1Q " 42 ( (M a /TeV)0 c - 2 ' < U » 
Comparing this with the estimate ([3]) one finds that the supernova physics is modified if 

A 2 

> 10" V 300 ms/At. (12) 



(M A /TeV) 2 

As was mentioned above, the best experimental bound on A 2 /M^ follows from the non-observation 
of oscillations of muonium atom [20]. We infer this bound from the review paper [6] taking into 
account that the definition of coupling constants in [6] and in the present work differ by the 
factor y/2: 

(iS7TwF <a2 - 90%CL (13 > 

From inequalities (|12p and ()13p it follows that even for At = 1 ms there is a room for a singlet 
bilepton which affects the supernova physics and evades direct experimental constraints. If one 
takes A = 0.1, then the bilepton mass should satisfy 

220 GeV <M A < 3.2 • (300 ms/At)" 1/4 TeV. (14) 



3 Note that in general the coupling matrix A = ||Ajj'|| is not invariant under rotations of the charged lepton 
basis. It transforms according to A' = U T AU = XU T U, where U is a matrix of a unitary rotation, U'U — 1. 
However, in case when U does not contain complex phases one gets = U T and A' = A. The role of complex 
phases is investigated in more detail in the next section. 
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4 Scalar triplet. See-saw type II. 



4.1 Cross sections 

Let us now discuss LFV mediated by a scalar triplet A. We are especially interested in the 
case when such triplet is responsible for the generation of Majorana neutrino masses through 
the see-saw type II mechanism [2]- BO- In this mechanism scalar triplet additionally couples 
to the SM Higgs bosons, this coupling producing a vacuum expectation value for the neutral 
component of the triplet. The neutrino masses are proportional to this vev. 

The see-saw type II Lagrangian contains two major ingredients, a scalar-lepton interaction, 

£ua = ^2 XwLfiT 2 ALi> + h.c, (15) 
l,V 

and a scalar potential, which in its minimal form reads 

V = —MftH^H + f(H^H) 2 + M|Tr(A f A) + ^-=(jiH T ir 2 ^ H + h.c). (16) 

V2 



Here 



Li = f j J 18 a doublet of left-handed leptons of flavor I = e, r, H is a Higgs doublet, jj, is 

a parameter with the dimension of mass. Note that the definition of coupling constant jl varies 
in different works by a numerical coefficient. Our definition coincides e.g. with that in [9]. 

Note that due to the anticommutation of the fermion fields 3x3 matrix A = ||A^'|| is 
symmetric, 

A T = A. (18) 
The vev of the neutral component of the triplet reads 

(A ) = (19) 
2V2MI 

where v = V2{H°) = 246 GeV. Due to the triplet vev neutrinos acquire the Majorana mass 
according to 

m = 2(A°)A, (20) 

where m = \ \mui\ \ is a neutrino mass matrix in a flavor basis. One gets that in the see-saw type 
II model neutrino mass matrix m is proportional to the coupling matrix A. Therefore now A can 
not be exactly diagonal as in the previous section, because m has non-zero non-diagonal entries 
(we know this since neutrinos oscillate). However, current experimental data on neutrino masses 
and mixing are incomplete and do not fix the ratio of diagonal and non-diagonal elements. In 
particular, it may appear that m is approximately diagonal, and thus it is possible to achieve 
the suppression of the processes with |AL/| = 1 compared to the processes with |AL;| = 2. This 
issue is considered in detail in the next subsection. 

Lagrangian (|15p gives rise to all four reactions ([T]), potentially relevant for supernova. Their 
cross sections in the center of mass reference frame read 



a(ee ^ w ) = (|A ee | 2 |A w | 2 /Ml)(l - ml/2E 2 e )^l - m\l 1 Ej E 2 e /2n, (21) 
a(ev e ^ ^) = (\\ ee \ 2 \\^\ 2 /Mi)(l-ml/4E 2 e ) 2 E 2 /4Tr, E e = E Ue , (22) 
a{v e v e -> = {\\ ee \ 2 \\^\ 2 /Mi)Elj2ir, (23) 
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Table 2: Constraints from low energy experiments. We extract numerical values from the 
review paper [6] taking into account that the definitions of coupling constants in [6] and in the 
present work differ by the factor \/2. 

o(y e v e v T v T ) = (\\ ee \ 2 \\ TT \ 2 /Mi)E 2 u j2ir, (24) 

where neutrinos and electrons are assumed to be extremely relativistic. 

Density and chemical potential of electrons are higher than those for electron neutrinos; 
also an additional factor 1/2 persists in eq. (|22|) . As a consequence, the first reaction in (pTJ) 
dominates over the second and the third ones provided the c.o.m. energy is considerably greater 
than the threshold 2rra M . The conditions on effective coupling constants, either of which ensures 
the appreciable rate of LFV in supernova core, now read 

|A ee A w | > 10 - 3v / 300 mg/At _J_ (25) 



M 2 - V / TeV 2> 

A fiR A 



, r , > 2.4 • 10 _3 V300 ms/At ? . (26) 

It is shown in the following subsection that Am„ ~ X TT for the interesting set of neutrino 
parameters. Therefore it is always easier to satisfy condition (j25[) than condition (|26p ; for this 
reason we concentrate on condition (|25p in what follows. 

If one assumes A ee ~ A MA1 = 0.1 and takes At = 300 ms, then, analogously to the case 
considered in the previous section, supernova appears to be able to test Ma < 3.2 TeV. To 
compare, the LHC can probe Ma < 1 TeV [7]. 

4.2 Neutrino mass matrix, rare decays and experimental constraints on TeV- 
scale see-saw type II 

LFV mediated by bileptons is constrained by non-observation of rare decays ([6|) and muonium 
conversion, see [5] [22] . The widths of rare decays and the rate of conversion are proportional 
to | \i x i 2 A; 3 ; 4 /M A | 2 , where l\ — I4 stand for the corresponding leptons. Experimental bounds are 
presented in Table 

Evidently the most stringent constraint stems from \i —> eee decay. Therefore in the present 
contribution we focus on this constraint. It may be shown that if this constraint does not 
contradict the condition (|25p . then other constraints also do not contradict it pp. 

The /i — > eee constraint does not contradict the condition (|25p if 



|A MM /A eM | = Im^/m^l > 40^/300 ms/At. (27) 

The neutrino mass matrix in the flavor basis, m, is obtained from the diagonal mass matrix 
through the transformation |25j 



m 



U* ■ diag(mi,m2,m3) • W, (28) 
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where U = \\Un\\ (I = e,/i,T, i = 1,2,3) is a PMNS neutrino mixing matrix: 

s 13 e- iS 




I I r,„ so,, I 1 1 

si 3 e iS 




x diag(e iai/ V a2/2 ,l), (29) 

Gjj — COS $ij j — S1U • 

Note that eq. (|28p implies m T = m in accordance with eqs. (|18p .(|20 p . The explicit expressions 
for the entries of m may be found in [7]; for the sake of completeness we provide this expressions 
in our notations]^ 

m ee = a c 2 3 + Si 3 m 3 e 2tS 

m M/J = mie" Ml (si 2 c 2 3 + si 3 e~ i,5 ci2S23) 2 + m 2 e~' ia2 (ci 2 c 2 3 - si3e~ J<5 si2.S23) 2 + m 3 cj 3 sl 3 
m TT = m 1 e~ iai (si 2 s 23 - si 3 e~ tS C12C23) 2 + m 2 e~ ia2 (ci 2 s 2 3 + s l3 e~ t& S12C23) 2 + m 3 c 2 Vi cl 3 
me^ = ci 3 [dsi 2 ci2C23 + si3e i5 S23("T-3 - ae" 2 ^)] 
"i er = ci 3 [-(isi2Ci2S23 + si3e* 5 C2 3 (m 3 - ae~ 2lS )\ 

= s 23 c 23 {-b + cf 3 m 3 ) - si 3 de" l5 si2Ci2(c| 3 - s| 3 ) + s 2 l3 ae~ 2l& s 23 c 23 

(30) 

Here we define the parameters with the dimension of mass 



a 


mie~ iai c\ 2 + m 2 e~ ia2 s 2 2 




b = 


mie~ iai s\ 2 + m 2 e~ ta2 c\ 2 


(31) 


d = 


m 2 e~ ia ' 2 — mie~ iai 





If the masses mi and m 2 are quasi-degenerate, m\ ~ m 2 , and Majorana phases a.\ and a 2 are 
approximately equal, ai ~ a 2 (this case is of particular interest, as is shown below), then this 
parameters reduce to 

a~6~e~ iQl mi (32) 

d^e- tai (m 2 - mi ) (33) 

There are nine parameters in the mass matrix: three neutrino masses rrii , three mixing angles 
9ij, the Dirac phase 5 and two Majorana phases a\ )2 . Experimentally measured (or constrained) 
quantities are [26] : 



#12^34°, 9 23 ~ 45° or 9 23 ~ 135°, % < 13°, (34) 

Amli ~ 0.8 • 10" 4 eV 2 , Am| 2 ^ ±2.4 • 10~ 3 eV 2 . (35) 
In addition, cosmological considerations |27| demonstrate that 

mi,m 2 ,m 3 < 0.35 eV. (36) 

Unknown quantities are phases 5 and ol\ 2, the mixing angle #13, the sign of Am 32 , the quadrant 
of #23, the absolute scale of neutrino masses. Depending on the sign of Am 32 two patterns of 
neutrino masses are possible: normal mass hierarchy (NH), m\ < m 2 < m 3 , and inverted mass 
hierarchy (IH), m 3 < m\ < m 2 . 

Formulas (|30|) allow to check whether condition (|27|) is valid for any specific choice of neutrino 
masses and mixing parameters. However, let us outline general trends of the dependence of 
\m^^/m efl \ on the neutrino masses, mixing angles and phases. 



4 It is worth noting that if the Majorana phase matrix would stay at the first place in the r.h.s. of eq. (|29|l . 
then the absolute values of entries of m would not depend on Majorana phases. The correctness of the ordering 
of matrices in eq. (|29f) is justified by the diagonalization of lepton mass matrices [25] . 
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Figure 2: Ratio \m^^/m efJi \ in case of exact CP-conservation (a± = ai = 5 = 0). If #13 = 
then |m w /m eA1 | does not depend on the quadrant of 623- Two thick horizontal lines correspond 
to |m MM /m eAt | = 40 and \m^^/m efl \ = 720, which are minimal values sufficient for LFV in 
supernova core to proceed at scales 300 ms and 1 ms accordingly. 

Generically 

\m^\ ~ max{m 2l m 3 }. (37) 

If one takes all three phases, ai,oi2,5, randomly from the interval [0, 2tt], then in the majority 
of cases 

\m efl \ ~ max{m 2 ,si3m 3 }, (38) 

which implies \m^^/m efJi \ < 50 Hence, in order to fit inequalities (|2"Tj) one should first of all 
impose some bounds on the phases. Note that if all three phases vanish (which correspond to 
the exact CP-conservation in the neutrino sector), m etl reduces to 

m e(U ^ si2Ci2C2 3 Am2i/(mi + m 2 ) + sis^Aml/ '(m 2 + m 3 ). (39) 

The first term here is small compared to rri2 if m\ ~ mi- The second term is small compared 
to 771,3 if either 7772 — 7773 or S13 is small. Therefore we come to a conclusion that small phases, 
large masses (which implies degeneracy of masses) and small #13 favour LFV in the supernova 
core. This is analogous to the conclusions of ref.[28], in which collider and low-energy LFV 
signatures were considered. 

We illustrate and refine this conclusions by means of plots depicted in Figs. On Fig. 
[2] the ratio (m^/me^l is plotted as a function of a minimal neutrino mass (mi for NH, 7773 for 
IH) for different values of #13 and different choices of the quadrant of #23 . All phases are taken 
to be zero, a\ = 02 = S = 0. One can see that in case of quasi-degeneracy of all three neutrino 
masses, m\ ~ 7772 — 7773 > 0.1 eV, condition (j27|) is fulfilled for any allowed value of #13. In case 
of a less stringent restriction, 777-1 — m 2 ^ 0.05 eV, condition (|27p with At = 300 ms is fulfilled 
for 0i3 < 3°. Note that two terms in m e ^ have different signs either in case of NH, C23 < or in 
case of IH, C23 > 0, see eq. ([39|) . This explains, particulary, the discontinuity at mi ~ 0.006 eV 
in case of normal hierarchy, #13 = 3°, 823 = 135°. 

To what extent is it possible to relax the requirement a\ = Q2 = 5 = 0? To answer this 
question we first note the following facts. 

• If 

at = a 2 , (40) 
To get this bound one should take into account that 7712 > 0.009 eV and in case of normal hierarchy 7713 > 0.05 

eV. 
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Figure 3: Scatter plots for the ratio \m^/m efl \. Each point corresponds to a random choice 
of neutrino parameters in the following ranges: m m \ a G [0.05 eV,0.35 eV],#i3 G [0°,3°],(5 G 
[-90°, 90°], a t G [-25 - 3°, -25 + 3°], a 2 G [-25 - 3°, -25 + 3°]. If the hierarchy is normal, then 
mmin = mi, otherwise m m ; n = 771,3 



then the parameter d reduces to e iai (7772 — ttii), as was pointed out above. If we for a 
moment put #13 = 0, then condition (|40p ensures that m e ^ is generically small: m etl = 
Cl3Cl2Sl2C23e" lai (7772 - 77li) Am^ /(mi + m 2 ). 

• If a more restrictive condition is superimposed, 

ai = a 2 = -25, (41) 

then m eil appears to be a sum of two generically small terms (with some relative phase 
factor) even for #13 7^ : 

m efl ~ ci 2 si2C23e 2i,5 A77i|i/(T77i + 777, 2 ) + s 13 e l5 s 2 3^ml 2 / ("72 + "2-3)- (42) 



Condition (|4ip is more general than the requirement a± = a 2 = 5 = 0, and embraces the 
latter as a specific case. We illustrate allowed deviations from this condition on FigJHJ One can 
see that equalities (j4"T|) should be valid up to ~ 3° for At = 300 ms in case when #13 < 3° and 
m min ^ 0.05 eV. One may further relax condition (|41|) imposing more restrictive requirements 
on the degeneracy of neutrino masses and smallness of #13 . On the other hand the specific value 
of oti ~ «2 — —25, the quadrant of 623 and the hierarchy of masses generically do not affect the 
value of 1 777^/777^1 significantly. 

Now it is easy see that in the interesting range of neutrino mass-mixing parameters the last 
two diagonal elements of the matrix A (and of the matrix 777) should be quasi-degenerate - a 
claim made in the end of section B~T1 This statement is verified numerically in pp. 

As was mentioned above, as far as the condition (|27p is fulfilled, not only /i — > eee constraint 
but all other constraints from the rare decays are satisfied. It is also worth mentioning that the 
t — > /jee rare decay is the most useful one among r — > III decays for the purpose of constraining 
the impact of LFV on the supernova physics pQ. The importance of the r —> [iee measurement 
for probing the see-saw type II model was emphasized in [29J . 



5 Possible consequences of LFV on supernova physics 

Before discussing the influence of LFV on supernova physics, we first give a sketchy qualitative 
picture of a standard collapse (cf., e.g. [30]) in section [5TT1 Then in subsection 15.21 we discuss 
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the collapse with LFV. The most intriguing effect of increase of neutrino luminosity due to LFV 
is discussed separately in subsection 15.31 



5.1 Standard collapse 

There are several stages in the standard core-collapse scenario where different neutrino processes 
dominate: 1) neutronization in the transparent core at infall before bounce; 2) bounce, followed 
by shock-breakout with the most powerful neutrino burst; 3) accretion onto growing proto- 
neutron star with a stalled shock and possible shock revival; 4) Kelvin-Helmholtz stage of 
cooling of a hot proto-neutron star. The latter stage is characterized by a slow contraction of 
the hot core on a timescale of many seconds with the release of gravitational energy through 
neutrino emission. 

Consider a polytropic equation of state 

P = Kp 1 (43) 

The fundamental critical value of 7 is 4/3 in Newtonian theory. If the average 7 < 4/3 over a 
stellar core, then this core is dynamically unstable. The only mass that may have a hydrostatic 
equilibrium for 7 = 4/3 in Newtonian gravity is the Chandrasekhar mass 

M Ch = 2.018(4tt) ( — ) . (44) 

This relation is derived in all classical textbooks on stellar structure. Physically, it is explained 
by a chain of simple arguments. The pressure P in general is the flux of momentum p, P ~ pnv, 
where n is the number density and v is the velocity of particles. For extremely relativistic (ER) 
particles P ~ pnc, and for degenerate electrons p ~ PFermi ~ 1/distance ~ n 1 / 3 . By definition 
of Y e - the electron number density per baryon, n = pY e /m n , where m n is nucleon mass. This 
implies P oc n 4 / 3 = (/9YJ=/m n ) 4//3 - which means 7 = 4/3 in relation (|43p . (The definition of 7 
is <91og P/dlog p\s-) 

We see that K oc (Y^/m n ) 4//3 in f|43[) . On the other hand, gravity requires the equilibrium 
pressure in the star center to be force/area oc GjyM 2 /i? 4 oc GjsfM 2 ' 3 p^ 3 . When equated with 
P = Kp 4 / 3 this fixes M Ch in 03]). Since K oc (y e /m n ) 4/3 one finds M Ch oc (Y e /m n ) 2 . When 
the coefficients are substituted accurately, one gets 

M Ch = 5.8y e 2 M . 

Curiously, since, in natural units Gn ~ m ~m > one ma Y obtain an elegant expression 

ml,Y? 



M Ch 



PV 



m 2 n 



where mpi is Planck mass. 

In general, instead of Y e there should be a sum over all particles producing pressure (e.g. 
e and v). And instead of the nucleon mass m n there may be a mass of another particle if the 
latter dominates in producing gravity. 

The contribution of baryons to the pressure is not significant until p > p n ~ 3 x 10 14 (nuclear 
matter density). The equation of state during collapse is almost polytropic with 7 slightly lower 
4/3 due to small corrections from nuclei. The value of K, which depends on Yl = Y e + Y v and 
entropy S, is nearly constant since neutrinos are trapped. 

Of course, at the bounce and later, in the center of the core the pressure is dominated 
by nuclear forces. However, the initial collapse begins when the pressure is dominated by 
electrons. The onset of the collapse in standard picture is determined either by the drop of 
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adiabatic exponent 7 below 4/3 due to iron dissociation or due to neutronization (or both, and 
in both cases the collapse begins due to the lowering pressure when the star core contracts). 

The dynamics of collapse in this case can be described by self-similar solutions [31] which 
show that a collapsing core separates into a homologous (v oc r ) central part, and an outer, 
supersonically infalling part. When the collapse is halted and bounces, a shock wave is produced 
at the sonic point. There is a self-similar solution with the shock as well 

A standard estimate for the free-fall collapse timescale [30] gives 



t H^) 1/2 " LW/2ms ' (45) 

where p\2 is the average density in units of 10 12 g/cm 3 . Below we use notation p^ for density 
in units of W k g/cm 3 . 

This estimate is not very realistic because it assumes zero pressure. In reality the force pro- 
duced by the pressure of ER degenerate leptons (~ n 4//3 i? 2 ) and the gravity force (~ GM 2 /R 2 ) 
scale according to the same law (oc R~ 2 ). Therefore the contribution of the pressure force 
remains finite and non-negligible compared to the gravity force during the entire collapse. A 
better estimate for the collapse timescale t c which agrees with numerical simulations can be 
derived from self-similar models [31} [32] : 

t c ~6t f . (46) 

The mass of the inner core turns out to be about 10% greater than the Chandrasekhar mass 
for the same value of K, and is thus proportional to the lepton fraction Y 2 . 

Another important timescale characterizes neutrino diffusion. Only when the diffusion is 
less rapid than the collapse will the lepton fraction be frozen for the duration of the collapse. 
The standard expression for the diffusion timescale is 

3R 2 Af 2 /y/ 3 <7 
td = ~ r ^ ~ • (47) 

If one takes the values corresponding to the inner core already after bounce, M ~ 0.6M & , p ~ 
P14, cr ~ 10~ 39 cm 2 , one gets 

t d ~ 10 s. (48) 

Kelvin-Helmholtz time is a bit larger than the diffusion timescale because td must be multiplied 
by a factor reflecting the total energy storage over the energy in neutrinos. 

Neutrino signal basically consists of three consecutive periods. The first one is the neu- 
tronization burst - a short (~ 10 ms) powerful outburst of neutrinos which emerges when the 
shock wave reaches the neutrino sphere. Neutronization burst takes place few ms after bounce, 
and contains almost exclusively electron neutrinos, which are created during neutronization of 
stellar matter. The second period coincides with the existence of stalled shock and lasts ~ 1 
s. The third period corresponds to the cooling of the proto-neutron star; it is characterized 
by an exponential decrease of neutrino luminosity and energy and lasts ~ 10 s. At the second 
and the third stages neutrinos and antineutrinos of all flavors are in approximate equilibrium 
at the neutrino sphere (not inside the core!) and therefore are emitted approximately in equal 
quantities. 

Neutrino oscillations take place above the neutrino sphere. In the layer ~ 100 km above 
the neutrino sphere collective oscillations occur (see e.g. |33j). In the outer stellar envelope 
more familiar MSW resonant flavor conversion take place (see e.g. |34] ) . Oscillations modify 
the flavor content of the neutrino signal. In particular, they may convert a fraction of electron 
neutrinos of neutronization burst into non-electron neutrinos. 
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5.2 Collapse with LFV 

The rate of LFV reactions ([1]) during the infall phase is negligible, with possible exception of 
the last few milliseconds. There are two reasons for this. First, the cross-sections (|21|M|24|) 
are small being proportional to the typical energies of e and v e (i.e., to chemical potentials of 
e and u e ). Second, the rates n e ^ Ue {av) are small when the concentrations n e and n Ve are low. 
Moreover, the first two reactions in ([1]) have thresholds proportional to muon mass m^. The 
reactions ([I]) may become important the few ms before bounce when the densities in the center 
approach the nuclear matter density. 

Let us estimate, for example, when ee —> pp process becomes important. The reaction is 
possible when E e > m«. At low density an electron may obtain such high energy only due to 
an exponential factor like exp(— m^/T), however T is always appreciably lower than m«. At 
higher densities electrons become degenerate, and we should estimate their chemical potential, 
taking into account that electrons are extremely relativistic: 

p e « ll(/9 10 y e ) 1/3 MeV. (49) 

We see that p e surpasses the muon mass at p ~ 3 x 10 13 g/cm 3 for Y e = 0.3 and the process 
ee —> pp sets in. However, if p e is only slightly above m^, only a small fraction of total electron 
density n e = pY e /m n contributes to n e {o~v) of this reaction: only electrons near and above the 
Fermi surface take part in the reaction. 

Nevertheless, the process of transformation ee — > pp starts decreasing the electron pressure 
already at pi% ~ 1, and becomes very important at P13 ~ 10, when the chemical potential of 
electrons is twice the muon mass and the bulk of electrons contributes to the reaction. At this 
stage we may safely use an estimate for the reaction rate 

-~n e M w 300pi4 V e s-\ (50) 

where we take E e = 2m^, A = 0.1, and Ma = 1 TeV (this gives u~2x 10~ 46 cm 2 ). Thus the 
timescale for the process ee — > pp is about 10 ms a few ms before the bounce. Our computations 
with the BOOM code show that the density p ~ 10 13 g/cm 3 is achieved in the center about 5 
to 7 ms before the bounce in good agreement with estimate (|46p . and there is 1 to 2 ms from 
p ~ 10 14 g/cm 3 until the bounce. 

The disappearance of electrons in ee — > pp reaction starts actively contributing to the 
decrease of electron pressure - hence to the dynamics of the collapse. Though the effect is 
rather small (in a couple of ms we lose only 10 to 20% of electrons) still all small deviations of 
7 from 4/3 may play a role. Moreover, the process of conversion will continue on later stages. 
It will lead to the decrease of electron (and electron neutrino) chemical potential and to the 
growth of entropy and temperature, since electrons are taken away well below the Fermi surface 
(cf. [351 M- 

The two-way LFV reaction ee ^ pp leads eventually to the establishment of equality p^ = 
p e for the chemical potentials. Muons become moderately degenerate since p^ — becomes 
somewhat larger than T. When the equilibrium is already established, the muon gas has 7 ~ 5/3 
(P ~ pn^v and v = p/m^, hence P ~ p 2 n^/m^ cx n^/ 3 /m^), which is harder than for ER 
leptons, although softer than for nuclear matter of high density. However, in the beginning of 
the transformation ee — > pp, when the reverse process is not yet effective, the pressure of the 
muon gas is much smaller than that of the electron gas due to the large mass (and therefore low 
momenta) of muons. Thus in the beginning of the ee — > pp conversion, which may start several 
ms before the standard bounce, a decrease of P e can not be compensated by the increase of P^, 
and we have an effect of reduction of total pressure similar to neutronization. This may lead to 
a faster collapse, hence to a stronger bounce. 

There is an opposite conclusion in ref . [37] . They claim that the reduction of Y e (accompa- 
nied by the reduction of P e ) leads to a weaker bounce. Note, however, that a one-zone collapse 
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model was used in [37] . It may be good for a uniform homologously contracting core where 7 is 

close to 4/3. Even when electrons are transformed into trapped electron neutrinos the gamma 
for lepton gas is close to 4/3 (all leptons are extremely relativistic). In this framework they give 
a standard estimate of the mass of the homologous core (HC), which is just approximated by a 
Chandrasekhar mass M Ch : M uc ^ M Ch a 5.8Y?M Q (a bit more accurate estimate following 
[3Tj [50] is Mhc — 1-1-^Ch)- However, if the core started collapsing and then Y e drops fast 
enough and non- uniformly (e.g., faster in the center, as in our case), then the picture may be 
more complicated. Effective 7 lower than 4/3 means faster collapse, hence a partial loss of 
homology. Quantitatively, the effect will be small if the process ee — > fifi is very slow. Anyway, 
here we see the qualitative difference with the standard case which may lead to the deformation 
of the homologous flow at infall and even to the formation of a second strong shock (because 
the innermost parts of the flow, those that are in the deepest potential well, have to accelerate 
due to the decrease of pressure). The modification of the shock wave naturally leads to the 
modification of the neutronization burst. We postpone the discussion of quantitative results for 
future research keeping in mind possible new dynamical phenomena. 

Another place where the reduction of 7 may be interesting is the shock after the bounce. 
There is a stage of more or less steady accretion through a stalled shock. When the pressure in 
the shock is dominated by leptons, and not by nuclear matter, the reduction of 7 means higher 
compression (the density jump in any strong shock is (7 + l)/(7 — 1). However, when the free 
path corresponding to LFV processes ([1]) is much larger than a hundred of km this effect will 
not be noticeable. 

The equilibrium between electrons and muons implies the equilibrium between electron and 
muon neutrinos which is established through charged current interactions. Also v e may be equi- 
librated with i/„ and v T directly through the LFV reactions ([1]). Anyhow, as a result the energy 
possessed by leptons is redistributed between two or three neutrino flavors, and degenerate seas 
of Vp, and probably v T emerge. The important consequences of such redistribution are discussed 
in the following subsection. 

5.3 Increase of neutrino luminosity due to LFV and consequences for the 
delayed explosion scenario 

How does the emergence of degenerate seas of non-electron neutrinos influence the neutrino 
transport? The key feature of and v T propagation in nuclear matter is that they do not 
take part in Charge Current (CC) interactions, in contrast with v e . The CC cross sections in 
the inner core are several times larger than Neutral Current (NC) cross sections. In the case of 
scattering on free nucleons 



In the dense supernova matter this cross sections are modified (see e.g. |38j). Moreover, while 
calculating mean free paths one has to take into account Pauli blocking for final particles, which 
is different for electrons and neutrinos. Nevertheless one expects that the opacity of the inner 
core for non-electron neutrinos is considerably lower than for electron neutrinos. This means 
that the neutrino luminosity (total emitted energy per second) increases in case of LFV. 

How does the increase of neutrino luminosity affect the explosion? Burrows&Goshy [39] 
and Murphy&Burrows |40] analyze this issue and conclude that such increase facilitates the 
explosion in the delayed explosion scenario. More specifically, they argue that for any given 
accretion rate M there exists a critical luminosity L Ue p e . If the actual luminosity exceeds this 
critical luminosity, the supernova explodes. This is illustrated by Fig. 1 in [39J and Fig. 17 
in |40] , In [39] the analysis is based on ID models, why in |40] 2D models are explored. The 
comparison shows that the critical luminosity for 2D models is 70% the critical luminosity for 
ID models. 

Thus one may expect that LFV in the inner core tends to facilitate the delayed explosion. 




(51) 
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To obtain a quantitative picture one has to make numerical simulations of the explosion in 
which v e comes to equilibrium with i/„ (and, possibly, v T ) in the inner core. However to get a 
rough idea of what will happen one may use a simple trick which does not require considerable 
modification of the supernova simulation code. Namely, one may use a code with a standard 
physics but with Charged Current interactions artificially switched off for neutrinos in the inner 
core. This will emulate the situation when non-electron neutrinos dominate the transport in 
the inner core. 

Remarkably, the relationship between the delayed explosion scenario and the decrease of 
neutrino opacities in the inner core was already studied by Burrows&Sawyer |41j . They do not 
consider any new physics. Instead, they take into account the density and spin correlations of 
nucleons and recalculate NC neutrino-nucleon scattering rates. They find that this scattering 
rates are suppressed by large factors around and above nuclear density, i.e. exactly in the in- 
ner supernova core. They conclude that supernova cores are more transparent than previously 
thought. Then they make some simplified protoneutron star cooling simulations with the ar- 
tificially decreased neutrino opacities. They find that the early luminosities are not altered, 
but the luminosities after ~ 500 ms are increased by factors that range from 10% to 100%. A 
brief but instructive discussion of this effect and of consequences for the delayed explosion is 
presented in Section VII of ref. [H] . 

To conclude this section, we suggest that equilibration between v e and v„ T in the inner 
core due to New, Lepton Flavor Violating Physics may increase neutrino luminosity and thus 
facilitate the delayed explosion. To quantitatively test this idea it is necessary to make a 
numerical simulation with new physics included which lasts at least one second after bounce. 
However, to get a rough idea of what happens one may merely switch off Charged Current 
interactions in the inner core without introducing any new physics in the simulation. 

6 Summary 

We have shown that the lepton flavor violation, generated by the exchange of TeV-scale scalar 
bileptons, may drastically modify conditions inside a core-collapse supernova. An important 
example of such bileptons is a scalar triplet responsible for the neutrino mass generation in the 
see-saw type II scenario. In this scenario the matrix of scalar-lepton couplings is proportional 
to the neutrino mass matrix. The latter is currently partially fixed by neutrino oscillations 
and cosmological data. Also the effective four-fermion couplings, generated by the triplet, are 
restricted by low-energy experiments. We analyze the body of data and demonstrate that there 
exists an allowed region of the see-saw II parameter space in which LFV is strong enough to lead 
to the thermal equilibrium between electron and non-electron species inside the collapsing core. 
Roughly speaking, LFV affects the supernova physics if the matrix of scalar-lepton couplings is 
approximately proportional to the unit matrix, A^/ — \5u>, and the effective coupling A 2 /M^ 
is greater than 10 -3 TeV~ 2 ~ 10~ 4 Gf. The former condition corresponds to neutrino masses 
greater than 0.05 eV and neutrino phases which satisfy the relations \ai + 2S\, |ct2 + 2<5| <C 7r. 

Equilibration between electron and non-electron species is a striking feature, absent in the 
standard picture of the collapse. It starts to settle not earlier than few milliseconds before 
bounce, when the density exceeds 10 13 g/cm 3 . Such equilibration lifts the chemical potentials 
of muons and non-electron neutrinos higher than 100 MeV, somewhat decreases the chemical 
potentials of electrons and electron neutrinos, leads to the increase of temperature and entropy. 
In addition to Fermi seas of extremely degenerate electrons and electron neutrinos, which are 
present in the standard picture, seas of extremely degenerate non-electron neutrinos and moder- 
ately degenerate muons emerge in case of LFV. This may affect the supernova physics in various 
ways. An elaborated study should be carried out in order to accurately and thoroughly describe 
all possible consequences. In particular, collapse simulations with new interactions included are 
necessary. In the present work we restrict ourselves by just outlining various possibilities. The 
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most intriguing one is the increase of the neutrino luminosity in the first second of the collapse 
due to the larger free path of non-electron neutrinos compared to electron neutrinos. Such 
increase is known to facilitate the explosion in the delayed explosion scenario. Moreover, it can 
be read out from the neutrino signal of a future galactic supernova. 

Other possible modifications of supernova dynamics are basically due to the decrease of 
pressure in the center of the core, where LFV reactions proceed effectively. Such decrease gives 
rise to a non-homologous stage of the collapse which may start few ms before the "standard" 
bounce. This may lead to modifications of the shock wave and even to additional shock waves. 
This, in turn, may lead to modifications of the neutronization burst and even to emergence of 
several neutronization bursts. 

We emphasize once more time that all effects listed above should be verified by profound 
supernova calculations. Once such calculations are done, it is possible to use the results in 
various ways. On the one hand, one can try to construct a successful supernova explosion model 
taking into account lepton flavor violation. On the other hand, one can extract parameters of 
the see-saw type II model and other models which include bileptons (or restrict such parameters) 
from the already persistent or future astrophysical data. 

As a finale remark we note that the discussed signatures of LFV in supernova are rather 
general and may be relevant in the context of various microscopic models, not necessarily those 
which include bileptons. In particular, the required LFV may be generated by the spin flips of 
Majorana neutrinos with large magnetic moment [42] or by the flavor changing neutral currents 
in SUSY with R-parity violation [43J. 
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